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Abstract 

We show that, in the heavy quark limit, the hadronic matrix elements that enter B meson decays into two light 
mesons can be computed from first principles, including 'non-factorizable' strong interaction corrections, and ex- 
pressed in terms of form factors and meson light-cone distribution amplitudes. The conventional factorization 
result follows in the limit when both power corrections in l/m& and radiative corrections in a a are neglected. We 
compute the order-a s corrections to the decays B c i — > ir + ir~ , Bd — > 7r°7r° and B + — > 7r + 7r° in the heavy quark limit 
and briefly discuss the phenomenological implications for the branching ratios, strong phases and CP violation. 
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We show that, in the heavy quark limit, the hadronic matrix elements that enter B meson decays 
into two light mesons can be computed from first principles, including 'non-factorizable' strong 
interaction corrections, and expressed in terms of form factors and meson light-cone distribution 
amplitudes. The conventional factorization result follows in the limit when both power corrections 
in 1/m.b and radiative corrections in a s are neglected. We compute the order-a s corrections to the 
decays Bd — * tt + ty~ , Bd — » 7r°7r° and B + — » ty + ty° in the heavy quark limit and briefly discuss the 
phenomenological implications for the branching ratios, strong phases and CP violation. 
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The detailed study of B meson decays is a key source 
of information for understanding CP violation and the 
physics of flavour. The interest in this field is rein- 
forced by the numerous upcoming experiments that will 
test crucial aspects of B decay properties with unprece- 
dented scope and precision. Among the large number of 
B decay channels, two-body non-leptonic modes, such as 
B — » 7T7T, B — > ttK etc., open a particularly rich field of 
phenomenological investigation. A theoretical treatment, 
however, is generally complicated owing to the non-trivial 
QCD dynamics related to the all-hadronic final state. 

In this Letter we describe important simplifications 
that occur in the limit m,b 3> Aqcd, when the b quark 
mass is large compared to the strong interaction scale 
Aqcd- We find that in this limit the hadronic matrix el- 
ements for, say, B — > t:tt can be represented in the form 



(7T7r|Q| J B)={7r|ii|5)(7r|i2|0)- 

1 + J2 r n< + 0(A QC D/m b ) 



(1) 



where Q is a local operator in the weak effective Hamil- 
tonian and j\,2 are bilinear quark currents. Neglecting 
power corrections in Aqcd and radiative corrections in 
a s , the original matrix element factorizes into a form 
factor times a decay constant (we call this conventional 
factorization). At higher order in a s this simple factor- 
ization is broken, but the corrections can be calculated 
systematically in terms of short-distance coefficients and 
meson light-cone distribution amplitudes. This is sim- 
ilar in spirit to the well-known framework of perturba- 
tive factorization for exclusive processes in QCD at large 
momentum transfer as applied, for example, to the 
electromagnetic form factor of the pion. An interesting 
consequence of (Q) is that strong interaction phases are 
formally of order a s or Aqcd/ttii, in the heavy quark 
limit. If this limit works well, the approach discussed here 
allows us to calculate these phases systematically; CP vi- 
olating weak phases can then be disentangled. Here we 
present a numerical analysis of B — > irir decay amplitudes 



based on the heavy quark limit. We also briefly discuss 
important power corrections, which should eventually be 
estimated in order to obtain a satisfactory phenomenol- 
ogy at realistic b quark masses. Details of the argument 
that leads to the factorization formula @ below will be 
explained in a forthcoming paper. 

The effective weak Hamiltonian describing B decays is 
given by § 



K 



eff 



Gp 

7f 



A p CiQi + C2Q2 + CiQ l 



i=3...6,8 



(2) 



where A p = V* d V pb . The Q t are local AB = 1, AS = 
operators, and Cj the corresponding short-distance Wil- 
son coefficients. We neglect electroweak penguin opera- 
tors and all terms not relevant to B —> irir decays. 

The essential theoretical problem for obtaining the 
B — * 7T7T amplitudes is the evaluation of the hadronic 
matrix elements (mr\Qi\B). Let 7Ti denote the pion that 
picks up the light spectator quark in the B meson, and 
7T2 the pion whose valence partons are supplied by the 
weak decay of the b quark. In the heavy quark limit 
both pions emerge with large energy mg/2 (in the B rest 
frame) . Power counting based on the asymptotic form of 
the leading-twist pion distribution amplitude shows that 
a leading-power contribution to the {it^Q^B) matrix el- 
ement requires both valence quarks of 7T2 to carry energy 
of order m b . The qq pair is ejected from the weak inter- 
action region as a small-size colour singlet object. As a 
consequence soft gluons with momentum of order Aqc d 
decouple at leading order in AQc/3/mf,, and ix^ can be 
represented by its leading-twist light-cone distribution 
amplitude. On the other hand, the spectator quark in 
the B meson carries momentum of order Aqcd and is 
transferred as a soft quark to 7Ti, unless it undergoes a 
hard interaction. The endpoint suppression of the pion 
wave function is not sufficient to ensure the dominance of 
hard interactions. [We adopt the point of view that for 
realistic b quark masses perturbative Sudakov suppres- 
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sion does not cut off soft contributions efficiently enough 
before one enters the non-perturbative regime.] There- 
fore 7Ti cannot always be represented by its light-cone 
distribution amplitude. At leading power in AQcu/mb, 
we find that the soft interactions can be absorbed into 
the B — > 7Ti form factor. On the other hand, any interac- 
tion of the spectator quark with the quarks of 7T2 is hard 
at leading power and can be written as a convolution of 
three light-cone distribution amplitudes. This discussion 
can be summarized by the factorization formula 

(HP>(<im\B(p)) = / B ^V) [ dxTi(x)^(x) 

Jo 

+ [ dSdxdyTf'fox^BWMx^M (3) 
Jo 

which is valid up to corrections of relative order 
h-QCD/rrib- Here f B ^ n {q 2 ) is a B — > ir form factor eval- 
uated at q 2 — m 2 w 0, and $„. (&b) are leading-twist 
light-cone distribution amplitudes of the pion (B meson), 
normalized to 1. The T i ' denote hard-scattering ker- 
nels, which are calculable in perturbation theory. Tf 
starts at 0(pr s )) at higher order in a s it contains 'non- 
factorizable' gluon exchange, including penguin topolo- 
gies, see the first two rows of Fig. [I] for the corrections 
at order a s . Hard, 'non-factorizable' interactions in- 
volving the spectator quark are part of T- (last row 
of Fig. |]). Annihilation topologies also exist, but are 
power-suppressed in AQCD/mb- The significance of the 
factorization formula is that all the non-perturbative ef- 
fects in the B — ► hit amplitudes can be absorbed into the 
form factor and the light-cone wave functions. 




FIG. 1. Order a s corrections to the hard scattering kernels 
Ti (first two rows) and T/ (last row). In the case of T/, the 
spectator quark does not participate in the hard interaction 
and is not drawn. The two lines directed upwards represent 
the two quarks that make up 7T2. 

The following comments are in order: 

(i) When a s corrections are neglected T/ 7 is zero and 
T£ is an x- independent constant. Conventional factor- 
ization in terms of the form factor and the pion decay 
constant is then recovered as a rigorous prediction in the 



infinite quark mass limit. The perturbative corrections 
are process-dependent, but calculable. Their inclusion 
cancels the scale-dependence of the leading-order factor- 
ization result. 

(ii) The infrared finiteness of the hard scattering am- 
plitude follows because the infrared divergences in the 
first four diagrams of Fig. [I] cancel in their sum. This 
cancellation is the technical manifestation of Bjorken's 
colour transparency argument |1. Colour transparency 
does not apply to hard gluon interactions. These, how- 
ever, are suppressed by a s and are calculable. 

(iii) The hard scattering contribution to the B — > it 
form factor is suppressed by one power of a s relative to 
the soft contribution, in which the B meson spectator 
undergoes no hard interaction. As a consequence the as- 
sumption that B — ► 7T7T can be treated entirely in the 
hard scattering picture of |IJ] would miss the leading con- 
tribution in the heavy quark limit. 

(iv) The decay amplitude acquires an imaginary part 
through the hard scattering kernels. In the heavy quark 
limit, the strong interaction phases can therefore be com- 
puted as expansions in a s . In terms of hadronic inter- 
mediate states that saturate the unitarity relation this 
implies systematic cancellations among many intermedi- 
ate states with potentially large individual rescattering 
phases. An estimate of rescattering effects on the basis 
of Regge theory is not compatible with the picture that 
emerges in the heavy quark limit. 

(v) The factorization formula (^) generalizes to the de- 
cays into a heavy-light final state, if the heavy particle 
absorbs the B meson spectator quark. In this case the 
second line in (|3|) is power-suppressed and only the form 
factor term survives. An expression of this form has been 
used by Politzer and Wise to compute the 1-loop correc- 
tions to the decay rate ratio T(B — > D*ir)/T(B — > Dir) 
[Q. The factorization formula does not hold for heavy- 
light final states, in which the light meson absorbs the B 
meson spectator quark, or for a heavy-heavy final state. 
In this case, conventional factorization can also not be 
justified. 

The result of an explicit calculation of the B — ► -kit 
decay amplitudes at order a s can be compactly expressed 
as (TTTT\n ef f\B) = G F /V2Y, P=u ,cK{™%\B), where 

T p = a^rnr) (ub)v-A ® (du)v-A 
+ a^i^n) (db)v-A <8> (uu)v-A 
+ a 3 (7T7r) (db)v-A <8 (qq)v-A 
+ al(irir) (qb) V -A ® (dq)v-A 
+ a 5 (7T7r) (db) V -A ® (qq)v+A 
+ al(-JT7T)(-2)(qb) S -p^{dq)s+p. (4) 

The symbol <£> is defined through (7T7r|ji (8 = 
(7r|ji|-B)(7r|j2|0). A summation over q — u,d is implied. 
Note that the term proportional to a^-nn) results in a 
power correction that should be dropped in the heavy 
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quark limit. We will comment further on this term be- 
low. 

Together with a\ (txtx) = a\(jfK) = and the leading- 
order coefficient ^(tttt) = Cq + C5/N, the QCD coeffi- 
cients a^(im) read at next-to-leading order (NLO) 



(5) 
(6) 
(7) 



<(-) = C 4 + -C3---^-C 1 + y C 3 
+ y (G 4 + C 6 )) In A + (^ (Sp ) _ Cx 
+ ^(0) + G ff (l) - - + G 3 
+ (3G W (0) + G 7r (s c ) + G^l)) (C 4 + C 6 ) + G,, 8 C 8 

^(tttt) = C 5 + lc 6 + G 6 (— i* 1 - 12). 



hadronic parameter As. Since 3>b(£) has support only 
for £ of order Aqcd/'ttib, is of order Kqcd- 

Writing the transition operator T p in terms of the 
QCD coefficients ^(nn) is a convenient notation for 
phenomenological applications. The notation general- 
izes the conventional parameters a\ 2 [PL which are seen 
to be process-dependent beyond leading order. We em- 
phasize that in the present context the ^(tttt) are not 
phenomenological parameters, but genuine predictions of 
QCD in the heavy quark limit. The Wilson coefficients Cj 
entering the of (rnr) are to be taken at NLO @], where we 
consistently drop terms of 0(a 2 ) in The physical 

amplitudes derived from (^) are independent of the renor- 
malization scale (fj,) and scheme through 0(a s ). The co- 
efficients ai(7r7r)-a5(7T7r) multiply scale and scheme inde- 
pendent matrix elements of (axial-)vector currents. Ac- 
cordingly for ai(7T7r)-a5(7r7r) the scale and scheme de- 
pendence in the Wilson coefficients Cj is canceled by the 
0{a s ) corrections in the hard-scattering amplitudes. In 
the case of ag(7T7r), a scale and scheme dependence re- 
mains, which is precisely the one needed to cancel the 
corresponding dependence in the matrix elements of the 

(9) 

(pseudo-)scalar currents, multiplying ag(7T7r) in (M). Be- 



(8) 



Here C F = (N 2 - 1)/(2N) and N = 3 (f = 5) is the 
number of colours (flavours). [Note that our definition 
of C\ and Ci differs from the convention of pj, where 
the labels 1 and 2 are interchanged.] The internal quark 
mass in penguin diagrams enters as s p , where s u = and 
s c = m 2 /m 2 . In addition we have used (x = 1 — x) 



(10) 



F = _121 n il-l8 + ^+/^, 
mb 

fi = [ dx g(x)$„(x), G v , s = [ dxG s {x)<$>„(x), (11) 
Jo Jo 

G„{s) = [ dxG(s,x)$ w (x), (12) 



with the hard-scattering functions 

1 — 2x 2 

g(x) — 3 lnrr — Sin, Gg(x) = — , 

1 — x x 



(13) 



G(s,x) = -4 duu(l-u)ln(s-u(l-u)x-ie). (14) 



The hard spectator scattering contribution is given by 



N U(0)m 2 B Jo 



\x^ 



(15) 



where f v (Jb) is the pion (B meson) decay constant, 
niB the B meson mass, /+(0) the B — > it form fac- 
tor at zero momentum transfer, and £ the light-cone 
momentum fraction of the spectator in the B meson. 
f^ 1 depends on the wave function $>b through the inte- 
gral L <2£3>b(£)/£ = itlb/Xb- This introduces one new 



sides the ln(/z/m;,) terms the hard-scattering amplitudes 
contain a scheme dependent constant, which we have ob- 
tained in the NDR scheme as defined in |J. This fixes 
the scheme to be used for the NLO coefficients Ci. 

At NLO the factorization coefficients ^(tttt) acquire 
complex phases, entering through the functions g{x) and 
G(s,x) in ( |l3| ) and (|l4|). Being of order a s , these phases 
are generically small, except in cases where the lowest 
order contribution is numerically suppressed. This hap- 
pens e.g. for ^(tttt). Physically, the phases arise from 
final state rescattering, which is due to hard gluon ex- 
change, and hence perturbative, in the heavy quark limit. 
The generation of strong interaction phases through the 
penguin function G(s,x) has been discussed many years 
ago |0] and is commonly referred to as the Bander- 
Silverman-Soni (BSS) mechanism. In the present ap- 
proach, the gluon virtuality k 2 — xm 2 B in the penguin 
diagram, which has usually been treated as a free phe- 
nomenological parameter, has a well-defined meaning. 
The ^-dependence of G(s, x) is convoluted with the pion 
wave function $ w (a;), leaving no ambiguity as to the value 
of k 2 . In addition we identify a further source of rescat- 
tering phases, represented by the function g{x). This ef- 
fect corresponds to hard gluon exchange between the two 
outgoing pions. Together with the BSS mechanism, it ac- 
counts for the complete asymptotic rescattering phases in 
B — ► 7T7T in the heavy quark limit. 

Another novel result is the existence of the contribu- 
tion from hard scattering involving the spectator quark 
in the B meson, expressed by in (^). This mech- 
anism is completely missed in phenomenological mod- 
els of factorization. It is particularly important for the 
small coefficient a^imr), where it leads to a sizable en- 
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hancement. Using f v = 131 MeV, f B = (180 ± 20) MeV, 
/+(0) = 0.275 ±0.025, A B = 0.3 GeV and the asymptotic 
pion wave function $ ff (x) = 6xx, we find f^ 1 w 6.4. The 
poor knowledge of the B meson parameter Xb makes this 
number rather uncertain. 

Numerical values for the a^(7T7r) are shown in Tabic Q, 
using the pole masses m& = 4.8 GeV, m c = 1.4 GeV, 
the MS masses fn t (fh t ) = 167 GeV, (m„ + fh d )(2 GeV) = 
9 MeV and Ai|L = 225 MeV as input parameters, a^{mr) 
multiplies the h-QCD I Tib-suppressed, but chirally en- 
hanced combination 



2m 2 



m&(/i)(m„(/z) + m d (fi)) 



1.18 [at fi = m b }. (16) 



In the following analysis, we give two results, one ne- 
glecting a^irir) as formally power-suppressed, the other 
keeping the leading-order expression for ag(nn). 

It is now straightforward to evaluate the decay am- 
plitudes and branching ratios. The latter are given by 
Br(B -> tttt) = T B /(W7Tm B ) • \A(B -> 7r7r)| 2 5', where 

5=1 for 7T7T = 7T + 7T~, 7T~7T and S =1/2 for 7T7T 

tb denotes the B meson lifetimes: t(B + ) 
r{B d ) = 1.56 ps. The decay amplitudes read 



7T 7T 

1.65 ps, 



A{B d -> vr+O 



Gp 1 

* "^f TO s/+(o)/7r|A c 



ttV) = *^|m|/ + (0)MA e 



• [R b e-^ «(7T7r) + oy(7T7r) + <(7T7r)r x ) 

-(aS(7r7r)+o§(7r7r)r x )], 

^(B d ^ 

• [i^e^ 7 (-o5(7nr) + ^(tttt) + a$(7rcr)r x ) 

-(aS(7r7T)+ag(7r7r)r x )], 

- tTtt ) = z^m|/ + (0)MA c | • 

• fl 6 /\/2 • e - * 7 (a?(7rer) + o^tttt)) . 



(17) 



(18) 



(19) 



Here i? 6 = (1 - A 2 /2)|K&/Kb|/A, where A = 0.22 is the 
sine of the Cabibbo angle, 7 is the phase of V* b , and we 
will use |F c6 | = 0.039 ± 0.002, \V ub /V cb \ = 0.085 ± 0.020. 
We find the branching fractions 



Br(B d -> TT+7T-) 

6.5 [6.1] • 10~ 6 
Br(B d -> 7r°7r°) = 

5.2 [7.7] • 10~ 8 
Br(B~ -> tt^tt ) 



r'T+0.09 [0.18] e 1 - 12 - 7 ' 6 - 7 ] 

r^ + 0.73 [1.11] e -<-i3r[i49]>= 
4.3 [4.3] • 10~ 6 , 



(20) 

(21) 
(22) 



where the default values correspond to neglecting a^nn) 
and the values in brackets use o^{-ktt) at leading order. 
While the predictions for the tt + it~ and ir~ir° final states 
are relatively robust, with errors on the order of ±30% 



due to the input parameters, the decay into ir°ir° depends 
very sensitively on the parameter Xb that controls the 
hard spectator scattering. If it is significantly smaller 
than 0.3 GeV, a branching fraction of order 10 -6 cannot 
be excluded. Eq. (^0|) can be converted into a result for 
the time-dependent CP asymmetry as a function of the 
CKM angle a. The direct CP asymmetry in the ir + ir~ 
mode is approximately 4% • sin 7. 

The approach discussed here allows us to formulate, for 
the first time, rigorous predictions of QCD for exclusive 
non-leptonic B decays in the heavy quark limit. On the 
other hand, as the dependence on the formally power- 
suppressed coefficient a^inr) demonstrates, the asymp- 
totic limit may be problematic at m b ss 5 GeV and the 
applicability of the theory has to be decided on a case- 
by-case basis. The most important power corrections are 
those that depend on the chirally enhanced combination 
(|ilf). The a s corrections to all such terms can in fact 
be identified. However, the factorization formula breaks 
down in this case, because the relevant twist-3 wave func- 
tions do not fall off fast enough at the endpoints. A 
detailed discussion of this point, and of its phenomeno- 
logical consequences, will be given elsewhere. 
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TABLE I. The QCD coefficients of (tttt) at NLO for three 
different renormalization scales \i. Leading order values are 
shown in parenthesis for comparison. 



: m b /2 



\x = 2m b 



CI2 (7T7T 

a\ (ttty 

as (ttty 
a p 6 (mv 



1.047 + 0.033i 

(1.038) 
0.061 - 0.106i 

(0.066) 
0.005 + 0.003i 

(0.004) 
-0.030 - 0.019i 

(-0.027) 
-0.038 - 0.009i 

(-0.027) 
-0.006 - 0.004i 
(-0.005) 

(-0.036) 



1.038 + 0.018i 

(1.020) 
0.082 - 0.080i 

(0.140) 
0.004 + 0.002i 

(0.002) 
-0.029 - 0.015i 

(-0.020) 
-0.034 - 0.008i 

(-0.020) 
-0.005 - 0.002i 
(-0.002) 

(-0.030) 



1.027 + 0.010i 

(1.010) 
0.108 - 0.064i 

(0.200) 
0.003 + O.OOlz 

(0.001) 
-0.026 - 0.013i 

(-0.014) 
-0.031 - 0.007i 

(-0.014) 
-0.003 - O.OOli 
(-0.001) 

(-0.024) 
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